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The absorption spectroscopy of calcium-doped argon clusters is described in terms of an effective 
diatomics molecule Ca— (Ar n ), in the framework of semiclassical vertical transitions. We show how, 
upon choosing a suitable reaction coordinate, the effective finite-temperature equilibrium properties 
can be obtained for the ground- and excited-surfaces from the potential of mean force (PMF). An 
extension of the recent multiple range random-walk method is used to calculate the PMF over 
continuous intervals of distances. The absorption spectra calculated using this single-coordinate 
description are found to be in good agreement with the spectra obtained from high-statistics Monte 
Carlo data, in various situations. For CaAri3 , we compare the performances of two different choices 
of the reaction coordinate. For CaAr37, the method is seen to be accurate enough to distinguish 
between different low-energy structures. Finally, the idea of casting the initial many-body problem 
into a single degree of freedom problem is tested on the spectroscopy of calcium in bulk solid argon. 



I. INTRODUCTION 

Polyatomic molecules can be electronically excited in 
a global or local way, depending on the chemical nature 
of the system, and on details of the excitation mecha- 
nisms, such as the characteristics or of the laser pulse. 
In many systems, the chemical bonds at various places 
of the same molecule favor local excitations. This sit- 
uation typically occurs in chromophore-doped rare-gas 
clusters, where photons in the visible range may excite 
the chromophore but not the rare-gas atoms. This fea- 
ture is especially useful because it allows a precise anal- 
ysis by comparison with the isolated chromophore spec- 
troscopy. Actually, even though the chromophore is lo- 
cally excited, the geometric distortions and the vibra- 
tions of the solvent atoms act as a perturbation on its 
spectroscopic properties. From a theoretical point of 
view, these many-body cffcc£s-|are usually treated either 
with simple approximations,!™ or convenientbf.by simu- 
lation. In the previous papers of this series,™ we have 
presented an alternative approach to simulation, where 
the Gaussian theory of absorption by Wadi and Pollakcl 
was combined with the superposition approximationEl in 
the quantum regime at finite temperature. This method 
was applied to CaAr„ clusters of various sizes n in the 
range 2 < n < 146. 

In these molecular systems, the calcium atom is en- 
ergetically more stable when located on the surface of 
the argon cluster, rather than in a fully solvated site. 
At moderate temperatures, below the cluster melting 
point, the calcium atom gets the ability to glide over 
the argon cluster. cl The heterogeneous cluster CaAr„ can 
then be thought as an effective diatomic molecule, where 
the argon cluster would be replaced by a large pseudo- 
atom. Similar ideas have been previously used in chemi- 
cal physics or condensed matter physics, as in the effec- 



tive interactions hetween colloids or between polymers, 
or in the GirifalcoO or Gspann-VollmaitO interaction po- 
tentials between clusters of carbon or rare-gas atoms, re- 
spectively. 

In this work, we investigate the possibility of treat- 
ing the photoabsorption spectroscopy of CaAr„ clusters 
in this effective diatomics picture. For this purpose, we 
need to characterize the interaction on the ground- and 
excited-state potential energy surfaces in terms of a sin- 
gle reaction coordinate separating the calcium atom from 
the argon cluster. The effective absorption spectrum can 
then be constructed from the effective potential curves. 
The finite-temperature properties of a many-body sys- 
tem are conveniently expressed as ensemble averages in 
the canonical ensembles. It is possible to treat the reac- 
tion coordinate separately, and to perform some averages 
on the remaining degrees of freedom. The next stage 
consists of averaging over the reaction coordinate itself 
using its suitable statistical weight. In this context, the 
statistical average on the many-body system is reduced 
into a potential of mean force (PMF) calculation that 
will appear as a Helmholtz free energy. Several methods 
le to get the computational solution to this 
Lincluding umbrella samplingt2l or constraint 
Lil or m o re recently the multiple-range random 
walk algorithm J 1 2 H 13 1 which we have used in the present 
work. 
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The paper is organized as follows. In the next section, 
we describe the general method, and the algorithms used 
to calculate the effective properties of the Ca-(Ar„) clus- 
ter as a function of an internal reaction coordinate. We 



then apply in Sec. HI the method to various cluster sizes, 
at several temperatures. In particular, the choice of the 
reaction coordinate may not be systematically obvious, 
and we show how more information can be obtained from 
a carefully chosen coordinate. For a suitable choice, ab- 
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sorption spectroscopy of calcium in bulk argon can also 
be studied, and some results will be presented in section 
[II. We finally summarize and conclude in Sec. IV. 



This defines a potential of mean force (PMF) VF(Co) ac- 
cording to: 



W (£o) - -k B T lnp(Co) 



(6) 



II. 



FREE-ENERGY PROFILES FOR GROUND- 
AND EXCITED-STATES 



Our system is a CaAr„ cluster, described using the 
Cartesian atomic coordinates R = {xi, yi, Zi}, where 
subscript will be used for the calcium atom. Wc 
first present the main ideas and approximations of the 
method. 

The photoabsorption spectrum of the cluster is calcu- 
lated in a semiclassical way by assuming vertical transi- 
tions (Condon approximation) between the ground state 
surface Vb(R) and several excited-state surfaces V4(R). 
At each configuration R, the absorption cross section 
CTfc(R) is proportional by the square of the transition 
dipole moment at this point, scaled by the transition 
energy. Excitation from the ground state surface ther- 
malized at temperature T leads to the unnormalized ab- 
sorption intensity given by the sum over excited 
states, I{lo) — ^2 k 1k(oj). The intensity of absorption 
from the ground state to state k is 



The absolute value of the PMF is meaningful only if the 
probability distribution p is normalized. It can be arbi- 
trarily shifted by any additive term Wo, provided that a 
factor e w °/ fcsT is incorporated in the calculated observ- 
ables. 

Building upon these definitions, we can introduce the 
partial average O(Co) of observable O by restriction on 
the coordinate £: 



d(£o) 



J Q(R)<S[g(R) - Co] exp[-T/ (R)/fc B r]rfR 
/ 5[£(R) - Co] exp[~V (K)/k B T}dR 



O can also be written 

/ 0(R)5[£(R) - Co] exp[-V/ (R)/fc s T]dR 



(7) 



O(Co) 



(8) 

where Zp(Co) plays the role of a restricted partition func- 
tion. By definitions of p and W, Eqn. (||) and (||), the 
global thermal average (A(cu, R)} over the whole configu- 
ration space is given by the one-dimensional average over 
the coordinate C: 



I k (w) = ±Js{hu- [V fe (R) - Vb(R)]}<7 fe (R)e- v °W/ fc * 



(A(oj,R)) = / A(C) 



,-W(t)/k B T 



( /C. 



(9) 



(1) 

where the partition function Z(T) is calculated on the 
ground state surface. For each excited surface k, and up 
to a constant factor, T^(lj) can be written as a canonical 
average on Vq: 



T fc (w) = {A k (w,R)), 



(2) 



with the notations 



A fc (u>,R) = a k (R)6{huj - [V fc (R) - Vb(R)]}; (3) 



and 



(O) = i J 0(R)e- y °( R )/ fc - T dR, 



(4) 



for any observable 0(H). In the following, and for clar- 
ity reasons, we shall drop the subscript k indicating the 
excited state. The single excited surface will be denoted 
V*. 

Let now assume that a reaction coordinate or an order 
parameter C(R) can be defined, which characterizes the 
overall location of ground state calcium with respect to 
the argon atoms. Possible expressions for £(R) will be 
discussed at the end of the present section. At a given 
temperature T, calcium has a probability p(Co) of resid- 
ing at coordinate Co given by the canonical average 



T dn, 

Importantly, Eq. (||) indicates that the PMF W(£) plays 
the role of a Helmholtz free energy term corresponding to 
the deformation of the system along the single coordinate 
C- We have thus mapped the initial many-body problem 
into a simpler one-dimensional problem, where the initial 
ground-state potential has been replaced by the PMF, 
and the instantaneous observable by its partial overage 
over all other, unrestricted degrees of freedom. 

In the present context, further assumptions are needed 
in order to get the full picture of an effective diatomics 
to describe the CaAr„ cluster. ^From expressions similar 
to Eq. (||) effective surfaces Vb(£) and V"*(£), as well as 
effective cross sections cf(C), can be constructed. The 
absorption intensity is then approximated as 

T(w) w I{uj) 



A further approximation can be done, which will be 
checked later. We replace the function exp(— W/ksT) 
by a Boltzmann weight over the effective ground-state 
surface Vq. Hence 

= C f 6{%uj - [V*(C) - V (0}}<r(Oe- m)/kBT <%, 



(11) 

where the proportionality constant C accounts for the 
normalization of Vq : 



p(£o) = (*[£(B.)-&]>. 



(5) 



C-i = j e-mtVtBT^ 



(12) 
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In the following, we will only consider normalized absorp- 
tion spectra, and C will be dropped. 

Using this extra approximation, Eq. (|ll|) thus obtained 
is the complete analogue of Eq. (0) for a single-coordinate 
system. In order to use it in a practical situation, we 
need to calculate all effective quantities over a continu- 
ous range of the reaction coordinate £, starting with the 
free-energy profiles Vo(£) and V'*(f). For this purpose 
we use the multiple range random walk method of Wang 
and LandauEj recently extended to the calculation of po- 
tentials of mean force and free-energies. IHi 

Briefly, we introduce a function g{£) initially set to 1 
in the range of accessible values of £, and we set s(£) = 
lng(£). A Monte Carlo simulation is carried out using 
the following Metropolis acceptanc&4jiile between the old 
Roid and new R ncw configurations :I13EJ 

acc(R old -» R now ) = min[l, exp(-AF/fc B T)], (13) 

with AF — F(R 11CW ) — F(R id) and the (Landau) free 
energy F(R) = Vb(R) - fc s Ts[£(R)]. After the new con- 
figuration Rnow is visited, s is updated: s[£(R nC w)] — > 
s[£(R ncw )] + In/, or equivalently, g is multiplied by /. 
Here / is a fixed quantity, initially set to 2-2.5. After 
a large number of MC steps, a new iteration m starts 
where a = In/ is reduced by taking a. rn +\ = a m /2. In 
this algorithm, the function T(£) = -fegTs(^) smoothly 
converges to the potential of mean force W(£)i up to an 
additive factor. 

Once it has been calculated, the PMF gives also access 
to the restricted averages A over a much wider range of 
£ than normally accessible. For this purpose each con- 
figuration in the Monte Carlo simulation is rewcighted 
by the factor exp{— Vy[£(R)]}. One should then get a 
uniform probability distribution for £. The flatness of 
this distribution can be used as a check of-jthe PMF com- 
puted using the Wang-Landau method.Ei The effective 
averaged observables A are then obtained with the usual 
reweighting formulas u 

We now turn to the choice of the reaction coordinate 
£(R) in concern with absorption spectroscopy. A suitable 
coordinate should have a clear geometrical meaning. A 
straightforward definition is the "local" distance d(R) 
between calcium and the closest argon atom: 



d(R) 



mm r; 

i>i 



(14) 



where fi denotes the position vector of atom i. This def- 
inition has the advantage that it reduces to the usual in- 
teratomic distance for the true diatomic molecule CaAr. 
Furthermore, this choice does not suffer any problem in 
a MC framework, in spite of the discontinuities that def- 
inition (|l4|) introduces in the Jacobian. In contrast, this 
would not be the case in a molecular dynamics simulation 
and an alternative better choice would be the distance 
between the chromophore and the center of mass (com) 
of the argon cluster: 



x(R) 



1 n 

n £ — ' 



(15) 



Let us notice that the latter reaction coordinate is not 
practical in several cases, such as calcium in bulk argon. 
However, it can sometimes bring extra information with 
respect to the choice of d above, Eq. (|l4]), as will appear 
in the next section. 



III. APPLICATIONS TO Ca@Ar 

We have applied the ideas developed in the previ- 
ous section to two CaAr n clusters, namely CaAri3 and 
CaAr 3 7. The 4s 2 ground state potential is modelled by 
simple pairwise terms, and the AsAp excited states are 
modelled using a Diatomic-In-Molecules (DIM) Hamilto- 
nian. All details about the potentials and the parame- 
terization are given in Ref. [lj. 

For all systems investigated here, the reference absorp- 
tion spectra were obtained from classical parallel temper- 
ing Monte Carlo, using 10 6 cycles after 2 x 10 5 initial cy- 
cles discarded for equilibration. For CaAr37, we did not 
use parallel tempering because only the local properties 
of the isomers were needed. The potential of mean force 
was computed from the multiple range random walk al- 
gorithm using 20 iterations of each 10 6 cycles following 
2 x 10 5 thermalization cycles. The parameter / was ini- 
tially taken as 2.5, and the number of bins in the his- 
tograms of the reaction coordinate was set to 1000. 

CaAri3 is an obvious choice for testing the above meth- 
ods, because the Ar 13 is icosahedral, hence nearly spher- 
ical, and because in this system calcium occupies a cap- 
ping site of the icosahedral cluster. This is a virtual sup- 
port to the effective diatomic picture. As seen in our 
previous papersJj'Q the calcium atom can jump into the 
icosahedral shell at moderate temperatures T ~ 30 K, 
but the spectroscopic signature of these isomerizations is 
rather weak. 

In Fig. |l| we have represented the effective potential 
energy curves (ground and excited states) of CaAri3 at 
T = 30 K, computed using the local distance d between 
calcium and argon defined by Eq. (|l4|) . We have superim- 
posed the (shifted) potential of mean force on the ground 
state effective surface. The two curves are rather close 
to each other, suggesting that approximating Eq. (|l0| ) by 
Eq. (|Il]) is correct. All curves have similar variations as 
in the Ca-Ar diatomics,Eil except a global shift to lower 
energies due to the additional argon-argon interaction 
energies. The excited state surfaces show a regular be- 
havior, but we also note the presence of a crossing near 
d ~ 10 ao which is also present in the Ca-Ar pair. The 
corresponding effective absorption cross sections (j(d) are 
displayed in Fig. |[ Except at small distances d < 6eto, 
they remain very close to their asymptotic atomic value. 
Hence the approximation (AB) A x B in Eq. ( |l0| ) 
should be satisfied. 

The photoabsorption spectrum calculated using the di- 
atomics picture and the effective potential surfaces is 
represented as a normalized histogram in Fig. |[ along 
with the result of classical Monte Carlo simulations. 
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FIG. 1: Effective potential surfaces of CaAri3 at T = 30 K, 
as function of the local distance coordinate d. Vo and V* 
are the ground- and excited-states surfaces, respectively. On 
the lower panel, the potential of mean force (PMF) is also 
represented, after shifting its asymptotic value to the ground- 
state potential. 



nificant signature of the isomers with calcium inside the 
icosahedral shell, the extra peak found in Ref. [|] being re- 
placed with a shoulder in the blue wing of the red peak. 
This prevents a clear identification of the possible isomers 
in the effective diatomics calculation. 
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FIG. 3: Normalized simulated (solid line) and effective, 
diatomics-like (histogram) absorption spectra of CaAri3 at 
T = 30 K, using the local coordinate d. 
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FIG. 2: Effective absorption intensities of CaAri3 at T ■ 
as function of the local coordinate d. 



30 K 



The agreement is remarkable, as far as the positions 
and widths of the peaks are concerned. The diatomics 
method introduces some additional noise due the unac- 
curacy in the calculated effective potential curves. This 
noise could be partly reduced by taking a smaller num- 
ber of bins in the interval of distances d. However, this 
would also require one to reduce the number of bins in 
the histogram of the absorption spectrum. 

As compared to our previous study using the super- 
position approximation^ the absorption spectrum calcu- 
lated using classical Monte Carlo does not show a sig- 



The CaAr37 duster has a large number of stable low- 
energy minima,S characterized by very different geome- 
tries. Its global minimum is decahedral, and the next 
most stable isomer is a Mackay-type icosahedron. Due 
to their different spectroscopic signatures at low tempera- 
ture, these two isomers provide a way of confronting the 
diatomics method to more detailed spectroscopic data. 
In Fig. |^ we have represented the effective ground state 
potential curves as a function of the local coordinate 
d, for the two decahedral and icosahedral isomers, at 
T = 5 K. At such low temperature, and using simple 
Monte Carkuthe cluster is expected to be trapped in its 
initial basinJa As can be seen from Fig. |, the icosahe- 
dral structure remains more stable than the decahedral 
minimum. The equilibrium distance is slightly shifted be- 
tween the two isomers, which results from the different 
local arrangements of argon atoms near the calcium site. 
Neither the effective excited state potential curves nor 
the effective absorption cross sections are shown in this 
figure, as they display the same variations as in CaAri3 
with the same reaction coordinate. The normalized ab- 
sorption spectra calculated with the effective diatomics 
approach are compared in Fig. to the reference Monte 
Carlo data. Besides the presence of residual noise, we 
find again a good agreement in the positions and widths 
of the absorption peaks for the two isomers. At temper- 
atures higher than 5 K, multiple isomers become popu- 
lated at thermal equilibrium^ but they are separated by 
large energy barriers, making the sampling difficult even 
within the Wang-Landau approach. In an experiment, 
such a situation would correspond to the presence of sev- 
eral stable isomers in the cluster beam. Turning back 
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FIG. 4: Effective ground state potentials of CaAr37 at T = 
5 K as function of the local coordinate d, for the decahedral 
global minimum (solid line) and the lowest Mackay icosahe- 
dral minimum (dashed line). 
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FIG. 5: Normalized simulated (solid line) and effective, 
diatomics-like (histogram) absorption spectra of CaAr37 at 
T = 5 K, using the local coordinate d. Upper panel: decahe- 
dral minimum; lower panel: icosahedral minimum. 



to calculations, rather than trying algorithmic alteration 
to treat all the isomers in one calculation, it seems more 
appropriate to perform calculations for each isomer and 
then average calculations with proper weights, in a spirit 
similar to the superposition approximation^ 

When the barrier between isomers is not as high as in 
these pathological systems, isomerization is in principle 
included in the calculated effective properties. However, 
the choice of the local coordinate d(R) may not be ap- 
propriate to reveal features associated with the presence 
of several isomers. A more suitable reaction coordinate 
would show distinct values depending on the isomers we 
want to separate. There is indeed some arbitrariness in 



the choice of d(R), which is very much guided by the 
kind of process we want to focus on. In CaAri3, the dis- 
tance d com (R) between the calcium atom and the center 
of mass of the argon atoms has essentially two values 
depending on the calcium atom being in a capping loca- 
tion over the argon icosahedron or inside the icosahedral 
shell. In the latter case, d com can actually take 4 differ- 
ent values, one for each icosahedral site, but the 4 values 
are close to each other (around 10 oo), and relatively far 
from the 13.5 ao value in the global minimum. 

The effective potential curves for CaAri3 at T = 30 K 
using the d com coordinate are represented in Fig. [|. We 
could not extend the sampling of this coordinate below 
8 ao or beyond 18 ao because of extensive noise. The ef- 
fective ground state potential clearly shows the presence 
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FIG. 6: Effective potential surfaces of CaAri3 at T = 30 K, as 
function of the center-of-mass coordinate d C om- Vo and V* are 
the ground- and excited-states surfaces, respectively. On the 
lower panel, the vertical dashed line defines the two different 
structures of the cluster, where the capping site is occupied 
either by calcium (A) or by argon (B). 

of two sets of minima, which can be attributed to the 
calcium (A) or argon (B) capping atom sketched in the 
lower panel of Fig. ^, depending on the outcome of a 
quenching procedure. Regions (A) and (B) are defined 
by the distance d com being larger or smaller than 12.01ao, 
respectively. The effective excited states potential curves 
are also strongly influenced by this change in coordinate, 
and their variations reflect the two stable minima near 
10ao and 13.5ao, respectively. Therefore the d com re- 
action coordinate provides a structural order parameter 
that can distinguish between the different parts of the 
configuration space we are interested in. Being able to 
differentiate isomers (A) and (B) allows one to calculate 
the separate contributions of each region to the global 
absorption spectrum. We have reported in Fig. [7] the 
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FIG. 7: Normalized simulated (solid line) and effective, 
diatomics-like (histogram) absorption spectra of CaAri3 at 
T = 30 K, using the center-of-mass coordinate i om . On the 
upper panel, the normalized spectra obtained by considering 
only regions A or B of the potential curves are represented 
(see Fig.|). 



total absorption spectrum and its individual contribu- 
tions from (A) or (B). The overall agreement between the 
global spectrum and the Monte Carlo data is again good, 
and we notice that the blue shoulder near 23900 cm -1 is 
larger than when using the local coordinate. This shows 
that the sampling of isomers (B) is more efficient, maybe 
even slightly too efficient with respect to the actual er- 
godic result. Longer simulations and a better sampling 
of region (A) would be needed to reduce this error. 

By carrying the integral ( |Tl] ) over each region (A) or 
(B) separetely allows for a spectroscopic distinction be- 
tween the isomers. It can be seen that region (B) is 
responsible for the blue wing at 23900 cm -1 . The global 
spectrum is the weighted sum of the two separate con- 
tributions from (A) or (B), and the statistical weight of 
(B) was found to be around 25%, in agreement with the 
study in Ref. || It is interesting to discuss further the 
shape of the spectrum for each isomer. We discussed in 
Ref. [l4| that the blue and red bands of spectra, similar 
to that shown in Fig. |?| for isomer A, are associated to 
exciting the 4p orbital perpendicular or parallel to the 
cluster surface, respectively. In the case of isomer A, be- 
cause of the outer location of calcium the two possible 
parallel alignments of the 4p orbital are almost degen- 
erated in average and lead to two merging bands. In 
contrast, calcium is closer to argon atoms in isomer B. 
As seen in Fig. || this corresponds to a larger the split- 
ting between the two lower (V*) potentials curves in the 
region that is accessible from the ground state isomer B. 
As a result, the red band in Fig. ^ is splitted into two 
components for isomer B. A splitted red band has actu- 
ally been observed experimentally in a slightly different 



context, the Bal6s 2 1 Sq — * GsGp^^Pi) excitation on large 
argon clusters Jl3 

We just have seen that the spectroscopic properties 
of the chromophorc are substantially affected by partial 
solvation in the argon cluster, that they appear as sensi- 
tive probe of the local environment of calcium and that 
they are satisfactorily described by the present effective 
diatomic picture. Going further in that direction, we 
consider a fully solvated Ca atom in an argon matrix 
as examined in the Monte Carlo simulations of Ref. 
The present system is now a calcium atom surrounded by 
107 argon atoms in a face-centered cubic lattice, at con- 
stant density p = f.562 x fO 3 kg.m -3 and temperature 
T = 20 K. Bulk argon solvent is simulated using periodic 
boundaries conditions in the minimum image convention. 
The reaction coordinate is again the local distance from 
calcium to all other argon atoms, and the simulations 
were carried out with the same statistics as for clusters. 
The effective ground- and excited-states potential curves 
are displayed in Fig. ||. Since calcium remains fully sol- 
vated in an equilibrium position, the ground-state sur- 

Ca@solid Ar 
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FIG. 8: Effective potential surfaces of Ca in solid argon at 
density p = 1.562 x 10 3 kg.m -3 and temperature T = 20 K, 
as function of the local distance coordinate d. Vo and V* are 
the ground- and excited-states surfaces, respectively. 

face roughly has a distorted parabola shape, and cannot 
tend to zero at large distances d. As we expect, the 
excited-states surfaces are fully degenerated, giviag rise 
to a dynamical Jahn- Teller effect in absorption.tJ The 
photoabsorption spectrum calculated from the effective 
diatomics approach is compared to Monte Carlo results 
in Fig. |i| The three excited potential surfaces lead to 
three absorption peaks, but they are quite broad, and 
the precise identification is less easy than with the MC 
data. Still we observe a surprisingly good agreement for 
the locations and overall widths of the peaks, suggest- 
ing that the effective diatomics method is adequate for 
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treating chromophore-doped inert systems with general 
shapes. 
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FIG. 9: Normalized simulated (solid line) and effective, 
diatomics-like (histogram) absorption spectra of Ca in solid 



argon at density p — 1.562 x 10' 
T — 20 K, using the local coordinate d 



m and temperature 



IV. DISCUSSION AND CONCLUSION 

Spectroscopy of polyatomic molecules with many de- 
grees of freedom is a difficult theoretical problem. Be- 
yond simple harmonic approximations^ the most con- 
venient way to describe correctly the absorption inten- 
sity relies on numerical simulations based on relevant 
Hamiltonians for the potential energy surfaces. However 
the complex, multidimensional character of the PES is a 
burden for simple physical interpretations of the results 
obtained by conventional methods. In this paper, we 
have presented a simple alternative approach based on 
the separate treatment of a single reaction coordinate, 
the many-body nature of the problem being thermally 
averaged into effective (free-)energy potential curves. By 
reducing the many-body problem into a single coordinate 
problem, the system is considered as a pseudo diatomics 
molecule, where the interactions implicitely depend on 
the thermodynamical conditions (temperature, but also 
pressure, and possibly chemical potential) and include 
the fluctuations within the averaged pseudo-atom. This 
point of view is best suited to the case of CaAr„ clusters, 
where the visible photoexcitation is localized on the sin- 
gle calcium atom. In heterogeneous clusters, extra reac- 
tion coordinates could naturally appear. For instance, in 
CaKrAr„, the two distances between calcium and kryp- 
ton, and between calcium and the argon cluster could 
be treated on a same footing. In such cases, the present 
investigation could be easily generalized to provide effec- 
tive hypersurfaces depending on several coordinates. In 



this picture the example chosen above would be consid- 
ered as a pseudo triatomics Ca-Kr-Ar ra . Similar ideas 
could be applied to atom-diatom reactions solvated by 
an inert cluster. 

Computationally speaking, the present method in- 
volves the preliminary calculation of the potential of 
mean force, which incorporates the statistical averages 
over all but one coordinates. The PMF is then used in 
a reweighting simulation to obtain the effective interac- 
tions and absorption cross section in a broad range of 
distances. The overall numerical cost is therefore heav- 
ier than in a conventional simulation, by one order of 
magnitude at the very least. However, the calculation of 
the PMF can be achieved conveniently, independently of 
the reaction coordinate chosen, using the multiple range 
random walk algorithm of Wang and Landaujio-recently 
adapted to the problem of free-energy profiles £3 There- 
fore the interest of the present approach mainly resides 
in the extra interpretations it provides. For instance, 
it could be used to assess or quantify the local charac- 
ter of the excitation. It also lays some ground for further 
spectroscopic investigations. In particular, more detailed 
studies of the excited states effective surfaces could be un- 
dertaken. The influence of temperature or the number of 
inert atoms on these curves, the related possible conical 
intersections offer examples of future research directions. 

On a practical point, we have largely used a local dis- 
tance coordinate d, namely the geometric distance (i.e. 
the lowest distance) between calcium and argon atoms. 
This order parameter is relevant to describe the CaAr„ 
cluster as a pseudo Ca-(Ar n ) diatomics, for which most 
of the methodology developed here was meant. How- 
ever, the separate effects of different isomers are essen- 
tially smoothed out and averaged, and it may be much 
more profitable to use other system-dependent coordi- 
nates to get a richer information about structure. In the 
case of CaAri3, we have shown that the distance between 
calcium and the argon cluster center of mass was more 
appropriate in this purpose. In other, more complex clus- 
ters, one should take profit of the known structural in- 
formation to design suitable reaction coordinates. For 
instance the bond order parameters Q 4, file, W4, or We, 
introduced by Steinhardt and coworkers, E!j could be used 
to distinguish between icosahedral, octahedral, or deca- 
hedral isomers. EZHlStJ One could then extend the present 
work to metallic clusters, where the energy levels can 
change upon global shape deformationJ23 and for which 
appropriate reaction coordinates could be the eccentrici- 
ties or the Hill- Wheeler parameters. 

Quantum vibrational effects are not accounted for in 
the present description, essentially because the corre- 
sponding treatment within path-integral finite temper- 
ature Monte Carlo is numerically demanding, and be- 
cause reference calculations are not available. However, 
at moderate temperatures, the present approach could 
be straightforwardly applied to quantum-corrected po- 
tential energy surfaces using the Feynman-Hibbs effective 
potentials .EH 
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At a more general level, the present approach gives the 
framework both to describe a reactive system in terms of 
deformations along a reaction coordinate and to calcu- 
late the corresponding energy variation. In that sense, 
the present calculation allows one to deal with the de- 
scription of a chemical reaction in terms of the transition 
state theory in situations where part of the reactive sys- 
tem itself acts as a thermal bath. In the present case the 
reactive system is a calcium atom plus an argon cluster. 



The latter plays the role of the thermal bath. The reac- 
tion is the simplest that might be considered: solvation 
of Ca by the argon cluster. The choice of reaction coor- 
dinate is not unique and a possible choice is the distance 
between Ca and the closest argon atom. Finally the rele- 
vant energy term when the system is free to evolve along 
the reaction coordinate is the potential of mean force that 
was define in the present work as a Hclmoltz free energy 
term. 
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